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ABSTRACT 

We present an X-ray analysis of the radial mass profile of the radio-quiet galaxy cluster A2589 
between 0.015 — 0.25 r v i r using an XMM-Newton observation. Except for a « 16 kpc shift of the X-ray 
center of the R = 45 — 60 kpc annulus, A2589 possesses a remarkably symmetrical X-ray image and 
is therefore an exceptional candidate for precision studies of its mass profile by applying hydrostatic 
equilibrium. The total gravitating matter profile is well described by the NFW model (fractional 
residuals < 10%) with c vir = 6.1 ± 0.3 and M vir = 3.3 ± 0.3 x 10 14 M Q {r mr = 1.74 ± 0.05 Mpc) in 
excellent agreement with ACDM. When the mass of the hot ICM is subtracted from the gravitating 
matter profile, the NFW model fitted to the resulting dark matter (DM) profile produces essentially 
the same result. However, if a component accounting for the stellar mass (M*) of the cD galaxy is 
included, then the NFW fit to the DM profile is substantially degraded in the central r ~ 50 kpc 
for reasonable M^/Ly- Modifying the NFW DM halo by adiabatic contraction arising from the 
early condensation of stellar baryons in the cD galaxy further degrades the fit. The fit is improved 
substantially with a Sersic-like model recently suggested by high resolution N-body simulations but 
with an inverse Sersic index, a ~ 0.5, a factor of ~ 3 higher than predicted. We argue that neither 
random turbulent motions nor magnetic fields can provide sufficient non-thermal pressure support 
to reconcile the XMM mass profile with adiabatic contraction of a CDM halo assuming reasonable 
Mx/Ly- Our results support the scenario where, at least for galaxy clusters, processes during halo 
formation counteract adiabatic contraction so that the total gravitating mass in the core approximately 
follows the NFW profile. 

Subject headings: X-rays: galaxies: clusters — dark matter — galaxies: clusters: individual (A2589) 



1. INTRODUCTION 

The properties of dark matter (DM) halos are a pow- 
erful discriminator between different cosmological mod- 
els. Of particular importance is the distribution of halo 
concentration (c v i r ) with virial mass (M v ; r ). For the con- 
cordance cold dark matter (ACDM) cosmology the mean 
c v i r varies slowly over a factor of 100 in M v j r , whereas the 
scatt er remains very nea rly constant fe.g.. lBullock et all 
120011 IKuhlen et al.ll2005l) . At fixed halo mass, the distri- 
bution of concentrations is expected to vary significantly 
as a function of the cosmological parameters, including 
an, n, and w, the dark energy equ ation of state (e.g. 
Pol ag et a,ljl200l IKuhlen et aWm^ . 

The radial density profiles of CDM halos are fairly 
well described between approximately 0.01 — 1 r v ; r (where 
r v ; r is the virial r adius) by the 2 -para meter "NFW" 
model suggested bv iNavarro et all <|1997t). More recent 
numerical simulations with higher resolution show that 
CDM halos deviate slightly from this average NFW pro- 
file because t he density slope cha n ges continuously with 
radius (e.g.. INavarro et all 120041: iDiemand et all 120041: 
iGraham &: et alJl2005[) . However, there is not yet a gen- 
eral consensus on whether the central density slope is 
shallower or steeper than the NFW profile. 
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The central regions of DM halos also reveal vital infor- 
mation about the interaction between the stellar baryons 
and the DM during halo formation. If the stellar baryons 
condense much earlier than the dissipationless DM, it 
is expected that the baryons will adiabatically com- 
press the DM halo away from its pure NFW form (e.g., 
iBlumenthal et alJl!986l iGnedin et al.ll2004|) . It has also 
been argued that heating of the dark matter by dynam- 
ical friction with member galaxies counteracts adiabatic 
compression and leads to a total gravitati ng mass profile 
consi stent wit h the pure NFW profil e (e.g.. lEl-Zant et alJ 
US see alsolLoeh fc Pe7hl^l2003l) . 

In theory, galaxy clusters are excellent sites to study 
DM profiles, particularly in their cores, because they 
are DM-dominated deep do wn to a small fraction of the 
virial radius(~ 0.007r v ; r ; e.g.. lDubinskill998l:lLewis et alJ 
120031) and several powerful techniques exist to probe clus- 
ter mass profiles. In practice, however, it has proven 
quite difficult to obtain precise, reliable measurements of 
the DM profiles in cluster cores. Stellar dynamical stud- 
ies suffer from the long-standing and well-known problem 
of velocity dispersion aniso tropv, and despit e occasional 
claims to the contrary fe.g. JSand et alJ2004|) . cannot ob- 
tain precise constraints without restrictive assumptions 
on the form of the velocity dispersion tensor, and there- 
fore of the mass profile. Since the weak lensing approx- 
imation breaks down within the central w 100 kpc of 
clusters, only those clusters with giant arcs produced by 
strong lensing can be used to map the core mass dis- 
tribution. Unfortunately, giant arcs are preferentially 
produced in clusters with significant substructure in the 
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core, making them much less suitable for comparison to 
the average relaxed cluster profile predicted by simula- 
tions. 

X-ray observations of the hot intracluster medium 
(ICM) are a powerful probe of the DM in galaxy cluster s 
(e.g., see recent reviews by Bu ot^l2004l lArnaud 12005'). 
Since the pressure tensor of the hot ICM is isotropic, 
and the gas traces the three-dimensional cluster poten- 
tial well within the virial radius, X-ray observations are 
a vital complement to gravitational lensing techniques 
that probe only the projected potential. However, the 
X-ray method requires the ICM to be in approximate 
hydrostatic equilibrium. Morphological studies of X-ray 
images have shown that a large fraction of nearby clus- 
ters (z < 0.2) do have regular image morphologies and 
appear to b e nea rly relaxed on 0.5-1 Mpc scales (e.g., 
[ Mohr et alJ 119951 iBuote fc Tsaj|ll996t Pones fe Formanl 
1999), though no cluster is observed to be (or is expected 
to be) devoid of evidence of disturbance. 

Precisely how departures from hydrostatic equilibrium 
translate to errors in mass estimates is as yet not well un- 
derstood in terms of quantifiable measures of the X-ray 
image morphology, projected temperature, or projected 
metallicity maps. What is clear is that errors in mass 
estimates should b e minimized for clusters without obvi- 
ous substruct ure l l Tsai et al.lll994| : |Bu^e^^^saJ ^95t 
lEvrard et all 1199(1 though see lHallman et alJ 120051. 
That is, clusters with regular X-ray isophotes, centrally 
peaked cores and no central AGN-induced disturbance. 
Further evidence for a relaxed hot ICM is suggested 
by a temperature profile that ris es with radius from a 
minim um in the cluster core fe.g.. lDe Grandi fc Molendil 
2002) and by metallicity profile that is peaked at the 
cente r and declines monotonically with increasing r adius 
fe.g.. lDe Grandi et al.ll200l iBohrineer et al.ll20o|) . 

Unfortunately, it has proven difficult to find cluster 
candidates that are relaxed both on large (« 0.5 — 1 Mpc) 
scales and within their cores because evolved, cool-core 
clusters typically have X-ray and radio disturbances 
withi n their cores believe d to arise from AGN feedback 
(e.g., iBirzan et alJ 120041) . A small number of radio- 
quiet clusters, previously classified as co oling flows, have 
been observed with Chandra; A 2589 l)Buote fc Lewisl 
l200l. A 644 iTBuote et a.U 12001 . A1650 and A2244 
ijDonahue et al.l|2005|) . None of these clusters displays 
the central minimum in the temperature profile charac- 
teristic of cool-core clusters or has any evidence for AGN 
feedback. The lack of AGN evidence for current feedback 
in these systems m ay arise from heating due to merging 
l|Buote et al.i r2005l or unu sually powerful AGN feedback 
events in the distant past ijDonahue et al.ll2005|) . 

The cluster A2589 (z = 0.0414) is observed to be one 
of the most morphologically regular clust ers in terms of 
its X- ray emission on a scale of ~ 0.5 Mpc ijBuote fc Tsail 
1996). The high-resolution Chandra image confirmed a 
highly symmetrical X-ray image morphology down to the 
center, although with some evidence for a small isophotal 
center shift of ~ 10 kpc. The highly regular X-ray image 
morphology combined with no evidence for AGN feed- 
back make this cluster an excellent candidate for studies 
of its mass profile by making the approximation that the 
hot gas is in hydrostatic equilibrium. Because of the low 
photon statistics of our shallow Chandra observation we 
were unable to place strong constraints on the mass pro- 



file, though an NFW profile was found to be consistent 
with the data within r ~ 150 kpc. 

In this paper we improve upon our previous study by 
using a higher quality XMM-Newton observation to ex- 
tend our analysis to ~ 400 kpc corresponding roughly to 
~ l/4r vir . The higher quality data will allow us to test 
different mass profiles proposed so far in the literature 
to account for the dark matter and the influence of the 
central bright dominant galaxy. 

Our assumed virial radius is defined as the radius of a 
sphere whose mean density is 104.7 times the critical den- 
sity of the universe. This value is estimated at the red- 
shift o f the cluster and for an (Cl m , fl\ ) = (0.3,0.7) uni- 
verse l)Brvan fc Normanlll998)) . (The redshift of A2589 
corresponds to an angular diameter distance of 171 Mpc 
and 1" = 0.82 kpc assuming f2 m = 0.3 and J7a = 0.7.) 

2. OBSERVATION AND DATA PREPARATION 

A2589 was observed by XMM-Newton for ~ 46 ks 
during revolution 822. The observation was performed 
in full window mode with the three EPIC cameras 
equipped with the thin filter. For the data reduction 
we used SAS 6.0, and for the spectral analysis we used 
XSPEC 11.3.1. 

The observation suffered from several periods of back- 
ground flaring. First, we excluded point sources that 
were identified by visual inspection. Then we isolated af- 
fected time intervals from visual inspection of light curves 
extracted in both high-energy (10-12 keV for the EPIC- 
MOS and 10-13 keV for the EPIC-PN) and low-energy 
(0.5-2.0 keV) bands. The light curves were extracted 
from regions far from the center of A2589. After excis- 
ing the flaring time intervals we arrived at cleaned expo- 
sures times of ~ 16.6 ks (MOS1), ~ 17.6 ks (MOS2), and 
~ 12.8 ks (PN). (Note only single events were used for 
the PN.) 

Since the cluster emission covers all the EPIC field 
of view it is not possible to obtain a local background 
estimate without accounting for source contamination. 
The standard procedure in this case is to use back- 
ground templates obtai ned from nominally blank fields 
ijRead fc Ponmanl 12003). But these standard templates 
are not appropriate in regions at large radii where the 
source emission is comparable to the background because 
of CXB spatial variations and, possibly, residual flaring 
particular to a given exposure. Consequently, we esti- 
mate the background by simultaneously fitting a model 
consisting of source and background components to spec- 
tra as far away from the source center as possible. 

Our procedure to subtract the backgroun d is an small 
updat e of the approach summar ized in |BijQ le_el_alJ 
1120041) that is explaine d fully in IBuote et alJ l|2006ft : 
iGastaldello et all 1)2006(1 . In sum, the model consists of 
components for (1) the cluster emission, (2) the Cos- 
mic X-ray Background (CXB) represented by two com- 
ponents for the soft thermal foreground emission and one 
power-law model for the extragalactic non-thermal back- 
ground, and (3) the quiescent instrumental background 
consisting of a continuous spectral component repre- 
sented by a broken power-law (whose parameters were 
determined from fitting the out-of-field of view events) 
and several gaussians to account for the instrumental flu- 
orescent lines. We model the non-thermal extragalactic 
CXB component due to unresolved AGNs using a power- 
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law with a slope T — 1.41. For galaxy clusters, because 
of their high gas temperatures, the contribution of the 
hot gas emission is partially degenerate with this high- 
energy power law component. Therefore, we required 
that the normalization of the CXB po wer-law component 
was in itially set to that obtained bv lDe Luca fc Molendil 
l)2004() a nd left free to vary only within the cosmic 
variance (jBarcons eTaIll20(i(y . All the cosmic compo- 
nents are absorbed by the Galactic column density value 
measured in the cluster r egion Njj = 4.15 x 10 20 cm~ 2 
( Dic kev fc Lockmanl lT990). Finally, we noticed that the 
broken power-law that we use to approximate the instru- 
mental background of the PN detector does not fit well 
the high-energy portions of the spectra in all of the an- 
nuli simultaneously. (See 21 for annuli definitions.) This 
may be ascribed to a residual flaring component of soft 
protons not completely removed that appears only in the 
inner r egions as it is focused an d vignetted by the satellite 
optics (Read fc Ponmanll2003fl . We were able to obtain 
an acceptable fit if we allowed the shape of the broken 
power-law in the outermost annulus to vary separately 
from the other annuli. 

Response matrices were generated for each region us- 
ing standard SAS tasks, which enable the generation of 
photon-weighted ancillary response (ARF) files. How- 
ever, these tools only enable redistribution matrix files 
(RMFs) to be generated appropriate for individual chips 
of the PN. Since we desire to analyze circular annuli that, 
in general, cover multiple chips, we created a photon- 
weighted RMF for each annulus combining the RMFs 
generated for appropriate ranges of CCD rows, between 
which the response is known to vary appreciably. This 
weighting was not performed for the MOS since the SAS- 
produced MOS RMFs do not vary over the field of view. 

3. IMAGING ANALYSIS 

We display the MOS1 image in Figure The image 
is not exposure corrected because the standard process- 
ing does not provided accurate exposure maps. Since 
the MOS2 and pn images have more chip gaps and 
bad columns, we only display the MOS1 image. The 
diffuse emission is remarkably regularly shaped, being 
moderately elliptical, and fills the entire field of view. 
We measure the centroid and the ellipticity of the X- 
ray surfac e brightness using t he m oment method de- 
scribed by iCarter fc Metcalfe ( 198Q() and implemented 
in o ur previous X-ray studies of galaxies and clusters 
('e.g. JBuote k, Canizareslir994j) . This iterative method is 
equivalent to computing the (two-dimensional) principal 
moments of inertia within an elliptical region. The el- 
lipticity is defined by the square root of the ratio of the 
principal moments, and the position angle is defined by 
the orientation of the larger principal moment. Follow- 
ing our previous s tudy of the elliptic ity of the Chandra 
data of NGC 720 (|Buote et al.l l2002t) we removed point 
sources and replaced them with smoothly distributed dif- 
fuse emission using the CIAO task dmfilth. Because 
this method for computing ellipticity cannot account for 
chip artifacts we restricted the analysis to the central 
chip on the MOS1. 

We defined the following annuli with outer semi-major 
axes (a) expressed in units of kpc (15, 30, 45, 60, 75, 
100, 150, 225). The ellipticity is rather uncertain within 



a ~ 30 kpc because of the small number of pixels (and 
resolution elements). For a > 60 kpc the ellipticity falls 
from near 0.30 to 0.25 at the outer radius with a position 
angle near 10° N-E over this entire range. 

The i ntensity weighted c enter shift, defined similarly to 
that in lMohr et al.1 1)1995|) . is 5.4 ± 0.5 kpc for these an- 
nuli. However, the maximum difference between annuli 
centers is observed between the annulus 45-60 kpc and 
the central circle (0-15 kpc) where we obtain a difference 
of 16 kpc in center positions with negligible statistical 
error. This 45-60 kpc annulus appears to be slightly off- 
set from all of the other annuli. Note that the X-ray 
peak (R.A.: 23h23h57.4m, DEC: 16d46m38s) coincides 
precisely with the optical position of the cD galaxy NGC 
7647. 

The most general non-rotating, self-gravitating equilib- 
rium configuration is the triaxial ellipsoid. The cquipo- 
tentials of a triaxial ellipsoid are concentric surfaces that 
themselves are nearly ellipsoids. Since in hydrostatic 
equilibrium the hot gas traces exactly the same shape 
as the gravitational potential regardless of the ra dial 
temperature profile (jBuote fc Canizares1ll994ll996|) . we 
expect a relaxed cluster to have concentric isophotes 
that are nearly elliptical in shape, though the elliptic- 
ity may vary with radius. (The position angle of the 
X-ray isophotes also need not be constant with radius.) 

The small center shift we measured above repre- 
sents a deviation from equilibrium, particularly near the 
isophotes with a = 45 — 60 kpc. We searched for higher 
order deviations from elliptical symmetry in the following 
manner. Using the ellipse package within IRAF-STSDAS 
we fitted elliptical isophotes to the MOS1 image; we re- 
stricted this analysis to the central CCD to avoid com- 
plications associated with the chip gaps. Then we con- 
structed a model image from the fitted isophotes and 
subtracted it from the raw image. The resulting residual 
image is also displayed in Figure Q We compared the 
model image and raw image in 4 sectors defined near the 
45-60 kpc annulus; i.e., inner radius 40 kpc, outer radius 
75 kpc. For each sector we obtain the following ratios 
of counts between the raw and model images: North 
1.01 ± 0.02, South 1.04 ± 0.02, East 0.92 ± 0.02, West 
0.94±0.02. These values represent quite reasonable scat- 
ter about unity considering uncorrected exposure varia- 
tions. We conclude the image does not show significant 
deviations from elliptical symmetry apart from the mod- 
est center shift noted above. 

We also searched for evidence of asymmetry in the tem- 
perature (and metallicity) structure using a hardness ra- 
tio map. Using the combined MOS1 and MOS2 images 
we constructed the map as an (S-H)/(S+H) image, where 
S and H are respectively images in soft (0.5-1.3 keV) and 
hard (1.3-8.0 keV) bands; the image was then smoothed 
with a gaussian kernel with width of 1 pixel. We display 
the hardness-ratio map in Fig. It is quite regular and 
shows no evidence for significant deviations from ellipti- 
cal symmetry. 

\ s first reported by iDavid et alJ <)1996l) there is a 
strong alignment between the optical isophotes of the 
cD, the X-ray isophotes, and the galaxy isopleths of 
the host (super-) cluster. It is also worth noting that 
the ave rage ellipticity of th e cluster galaxy isopleths is 
~ 0.3 llPlionis e t all 119 911). consistent with the X-ray 
isophotes. iBeers et al.l (1991) reported an offset of the cD 



4 



Zappacosta et al. 




Fig. 1. — {Left) MOS1 image in the 0.5-8.0 keV band displayed with logarithmic intensity scaling. This image has not been exposure 
corrected due to the inaccurate exposure maps available from the standard processing. No additional smoothing has been applied to 
this image. However, smoothed intensity contours are overlaid to guide the eye. (Right) Residual image of MOS1, displayed with linear 
intensity scaling, obtained by subtracting from the raw image a model constructed by fitting perfect elliptical isophotes (see text). The 
region displayed corresponds to the boxed area of the image on the left. The sectors indicate regions analyzed to assess non-ellipsoidal 
fluctuations as discussed in the text. 
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Fig. 2. — Hardness ratio image of A2589, displayed with loga- 
rithmic intensity scaling, generated using the 0.5-1.3 keV and 1.3- 
5.0 keV combined MOS1 and MOS2 images (see text for details). 
Contours represent the X-ray isophotes. 



velocity (~ 250 km s _1 ) from the other member galax- 
ies. However, these sparse data were deemed insufficient 
to obtain a robust result, and, consequently, this sys- 
tem wa s not i nclud ed in the systematic study of cD off- 
sets bvlBirdl II1994D. Ind eed, the kinematic al study by 
IGebhardt fc Beersl (fl991^ using the data of iBeers et alJ 



(1991) obtains only weak evidence for substructure. 

Evidently, A2589 is nearly in hydrostatic equilibrium. 
This is indicated by the strong similarity in the distribu- 
tions of galaxies and ICM, and because A2589 has one 
of the most regular X -ray images on the « 0.5 Mpc scale 
l)Buote fc Tsail ll996j) . The only significant asymmetry 
we detect is the small center offset associated with the 
45-60 kpc annulus mentioned above. We will consider 
whether this region displays anomalous behavior in our 
analysis of the azimuthally averaged spectral (and mass) 
properties below. 

4. SPECTRAL ANALYSIS 

Since the X-ray image and hardness-ratio map do not 
display substantial azimuthal variations, we focus our 
analysis on the X-ray spectral properties derived from 
circular annuli. We neglect the ellipticity of the X-ray 
isophotes since determining the flattening of the mass 
distribution is beyond the scope of this paper. (In any 
event, using elliptical annuli does not produce qualita- 
tively different profiles of the gas density, temperature, 
or metallicity when the results are expressed in terms of 
the mean radius R = a^fq, where a is the semi-major 
axis and q is the axial ratio.) 

We extracted spectra in the 0.5-8.0 keV energy band 
from all of the EPIC CCDs in a series of concentric circu- 
lar annuli centered on the X-ray peak. These annuli were 
originally constructed to have at least 8000 background- 
subtracted counts in the MOS1 and were defined to be 
larger than the XMM PSF and provide temperature con- 
straints of similar precision. Since the central radial bin 
has the highest S/N, we further divided it into two annuli 
(each still being larger than the XMM PSF). This allows 
us to probe the spectral properties down to ~ 0.015r v j r . 
The annuli definitions are listed in Table □ We do not 
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TABLE 1 

Spectral fit parameters for each annulus 



Annulus 


R 


R 


T 


^Fc 


^b 


^bi 


norm 


v 2 /dof 


Null Hyp. 




(kpc/arcmin) 


(kpc/arcmin) 


(keV) 


(solar) 


(solar) 


(solar) 


(10- 3 cm- 5 ) 




prob. (%) 


i 


0/0.00 


32/0.65 


3.26 ±0.08 


1.33 ±0.11 


0.55 ±0.25 


1.01 ±0.20 


1.38 ±0.04 


535.3/487 


6.39 


2 


32/0.65 


57/1.16 


3.45 ± 0.08 


0.95 ±0.07 


0.17 ±0.16 


0.47 ±0.19 


2.10 ±0.04 


654.6/606 


8.38 


3 


57/1.16 


93/1.89 


3.61 ±0.07 


0.78 ± 0.06 


< 0.22 


0.30 ±0.16 


3.42 ±0.08 


752.0/730 


27.8 


4 


93/1.89 


136/2.76 


3.40 ± 0.05 


0.77 ±0.06 


0.14 ±0.09 


0.16 ±0.08 


3.72 ±0.07 


772.9/716 


6.91 


5 


136/2.76 


193/3.92 


3.36 ±0.06 


0.54 ±0.05 


0.22 ±0.12 


0.52 ±0.09 


3.99 ±0.06 


653.6/688 


82.0 


6 


193/3.92 


289/5.87 


3.52 ±0.07 


0.57 ±0.07 


< 0.12 


< 0.14 


4.56 ±0.09 


726.4/676 


8.53 


7 


289/5.87 


489/9.93 


3.04 ±0.15 


0.43 ±0.07 


0.30 ±0.15 


0.28 ±0.19 


5.40 ±0.16 


493.5/450 


7.12 



Note. — Parameters derived by fitting an absorbed APEC plasma model to the spectra extracted in each annulus. Columns Zp e , Zq , Zsi show the 
abundance values obtained leaving free to vary Fe, S and Si. The "norm" parameter is the emission measure for the APEC code as defined in XSPEC. 
The last column report the Null Hypothesis Probability for the fit. The total 0.5-8.0 keV luminosity of these annuli is 1.47t° *6 1q43 cr S s_1 ( 90 % 
conf.). 



list results obtained for the outermost annulus (R > 9.93' 
out to the edge of the fields) that were used in the back- 
ground determination summarized in |21 At these large 
radii the source emission contributes little to the EPIC 
spectra, and thus small changes in our adopted back- 
ground model translate to relatively large changes in the 
derived source parameters. (We assess the systematic 
error resulting from the background level on our mass 
results in jJEI) We also restricted the upper energy limit 
to 5.0 keV in the penultimate annulus (i.e., the last an- 
nulus listed in Tabled because of excess hard emission 
that we could not remove with our background model. 

In each annulus we fitted a model consisting of an 
optically thin hot plasma (apec) modified by Galactic 
absorption. The free parameters are the normalization, 
temperature, iron, silicon, and sulfur abundances. All 
other elemental abundances are fixed so that their ra- 
tio with res pect to iron is solar assum ing the abundance 
standard of iGrevesse fc Sauvall (1998). Error estimates 
are obtained by simulating spectra based on our best- 
fitting models. From 20 Monte Carlo simulations we 
compute the standard deviation of each parameter which 
we report as the la error. (All quoted errors are la un- 
less stated otherwise.) 

The results are listed in Table [2 for the model param- 
eters, and we display the MOS1 spectra and best-fitting 
models for two annuli in Figure 01 The quantities we 
derive for the hot gas in each annulus are average values 
weighted by cluster emission projected along the line-of- 
sight. Below we will obtain the gas density and temper- 
ature as a function of three-dimensional radius by pro- 
jecting models along the line-of-sight and fitting them 
to the results obtained as a function of projected radius 
listed in Tabled We found this procedure to be appro- 
priate given the small number of radial bins and quality 
of our data. In [0 we also consider the effects on the 
derived mass profile of first deprojecting the data using 
the well-known onion-peeling method. 

We mention that the generally sub-solar ratios of Si/Fe 
and S/Fe abundances imply that most of the iron (~ 
80% -95%) in the hot ICM has been provided by Type la 
supernova, with convective deflagration explosion models 
favored over delayed-detonation models. These results 
are consistent with those we have obtained for galaxies 



and groups fe.g-. lBuote et all2003btlHumphrev fe Buotd 

5. GAS DENSITY PROFILE 

Our analysis of the gravitating mass and dark matter 
requires that we evaluate derivatives of the gas density 
and temperature with respect to the three-dimensional 
radius ( H7.1l and ij7.3|l . To reduce the noise in the deriva- 
tive calculations we fit simple analytic functions to the 
entire radial range of the gas density and temperature 
data. For the gas density (p g {r)) we integrate the quan- 
tity, p^A(T, Zpe), along the line-of-sight, where A(T, Zp e ) 
is just the APEC model with norm = 1. The temperature 
profile model is fixed to the best-fitting result obtained in 
the following section. (The iron abundance profile is also 
fitted with a simple model and extrapolated as a power- 
law outside the last data point. It also remains fixed 
during the fit.) The limits of integration are projected 
radius R and a maximum value, the latter of which we 
set near 2 Mpc corresponding to the inferred virial radius 
of the cluster below. (The results are not sensitive to this 
choice.) Consequently, the integral is obtained as a func- 
tion of projected radius R which we then evaluate over 
the radial width of each circular annulus defined in Table 
2] The result for each annulus is divided by A(T, Z-p e ), 
where now T and Zp e are the emission-weighted results 
quoted in Tabled for the annulus in question. 

As is standard for such studies, we begin by fitting 
the gas density profile with the well-known model 
<|Cavaliere fe Fusco-Femianolll978l). 



Pa = Pgo 



1 



(i) 



where p g o is the central gas density, r c the core radius, 
and the asymptotic slope is —3/3. The result of fitting 
this model to the data is shown in Figure 0] and the 
parameters are listed in Table |3 Although the single 
(3- model is not formally an acceptable fit, it provides a 
good representation of the radial profile. The fractional 
residuals are < 5% in all annuli except the center where 
a ~ 10% deviation is observed. 

Formally acceptable fits can be obtained by employing 
conventional modifications of the singl e model. First, 
we examine the "cusped (3 model" (jPratt fe Arnaudl 
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Fig. 3. — MOS1 spectra accumulated within (Left Panel) annulus #1 and (Right Panel) annulus #7. Each spectrum is fitted with an 
APEO plasma model modified by Galactic absorption as discussed in J2] Note that these models were fitted simultaneously with the MOS2 
and pn data, but we display only the MOS1 for clarity. 



TABLE 2 

Gas Density and Temperature Models 



Projected Density 
/3 model double /3 model cusp /3 model 


Projected Temperature 
Power-law model model 


/3 = 0.54±0.01 j3 = 0.57 ±0.02 /3 = 0.59±0.02 

r c = 76.5 ± 2.4 kpc r c i = 21 .9 ± 14.5 kpc r c = 119.2 ± 12.5 kpc 

p g o = 1.40 ± 0.03* r c2 = 95.2 ± 14.0 kpc ct c = 0.38 ±0.06 

p g0 ,i =2.1 ±8.0* p g ,c = 0.82 ±0.10* 

Pg0,2 = 1.1 ±0.1* 

X 2 /do/ = 27.6/4 X 2 M>/ = 0.67/2 x 2 /*>/ = 1.08/3 


Tioo =3.42 ± 0.03 keV Ai =3.6 ±1.5 
a p = 0.0±0.02 A 2 = 3.5±0.9 
r p = 100(fixed) 
ai =0.21 ±0.21 
ol2 = -0.09 ± 0.15 
7 = 0.45 (fixed) 
X 2 /dof = 21.3/5 X 2 /dof = 15.0/3 


Note. — Models of the gas density and temperature as a function of three-dimensional radius. The functions and 
parameters are defined in fj^jand The gas density normalizations (*) arc expressed in units of 10 — 26 g cm -3 . 



2002; ILewis et al.l l2003). 



P 9 = Pg,c 



2 3/3/2-a c /2 ( j_ 
Tr. 



1+1- 



I/3+- 



( 2 ) 

where the exponent a c is the slope of the power-law cusp 
at small radii and p g<c = p g (r c ). Although this model 
introduces only one additional free parameter over the 
single model, as is seen in Figure 0] and Table |2 the fit 
is excellent: residuals at all radii are reduced to < 1%. 

Second, we also investigate d adding a second model; 
i.e., a "double-/3" model (e.g. JXu et alJl998tlMohr et all 
1999). Even if we require the values of both compo- 
nents to be the same, the fit is excellent - just as good as 
the cusped model (see Figure0|and Table[2l . Since this 
model does not improve the fit over the that achieved by 
cusped model, but it introduces another free parame- 
ter, we will use the cusped model as the default profile 
in the mass analysis below. 

The values we obtained for are larger than we 
obtained previously for A2589 with Chandra data 
ijBuote & Lewisl 120041) . Using a single model, which 
was all that was required by the low-quality ACIS-S ob- 
servation, we obtained = 0.39 ±0.04; i.e., w 4<r smaller 
than the result obtained for the single model with the 
EPIC data. This difference is not surprising since the 
Chandra data were only fitted out to w 150 kpc com- 



pared to w 500 kpc in this paper. The values obtained 
from XMM agree well with the value of 0.57 quoted by 
IDavid et all <|1996|) obtained from ROSAT data. 



6. GAS TEMPERATURE PROFILE 

We display the temperature profile in Figure [S] corre- 
sponding to the values in Table ^ The approximately 
isothe rmal profile is s i milar to that measured by Chan- 
dra in lBuote fc Lewisl l)2004l) and does not s how the evi- 
dence for a cool core indicated by ROSAT ijDavid et alJ 
199(|. The lack of a cool core is unusual for a cluster 
having such a regular X-ray morphology, though the cen- 
trally peaked iron abundance pr ofile is similar to those 
found in cool core clu sters (e.g.. lDe Grandi et al.l 120041 
Bohringer et~alll2004l) . 

Starting with a model for T(r) we projected the quanti- 
ties TpgA(T, Z Fc ) and p 2 g A(T, Z Fe ) along the line of sight. 
(Here we use the best-fitting model for p g (r) obtained 
above, and the Zp c profile is fixed as before.) These 
quantities were evaluated within annuli in the same man- 
ner as done for the gas density. The projected emission- 
weighted temperature is obtained by dividing the first 
term by the second term. 

Since the temperature data are nearly isothermal, we 
initially fitted a single power-law, 



T(r) = Tj 



too 



( 100 kpc) 



(3) 
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Fig. 4. 



Radius (kpc) 
The radial profile of the projection of p\ ?. The data 



points arc proportional to norm/ A, where norm is listed in Table 
HI and A is the area of the annulus. The models have been fitted 
by weighting the projection integral by the plasma emissivity (see 
However, for display purposes we do not bin the models in the 
plot. Because the statistical error bars are very small, we enclose 
the data points with open circles. Each data point represents a 
radial bin as defined in Table HI The horizontal bars represent the 
width of the annuli. The solid line is the best fit cusp (3 model, 
the dotted red line is the single [3 model and the dashed blue line 
is the double (3 model. The cusp (3 model and double (3 model 
almost entirely overlap. The inset shows an enlarged view of the 
inner region. In the ratio (data/model) plot, red triangles, blue 
diamonds and black squares refer respectively to (3 model, double 
f3 model and cusped (3 model. 



where Tioo is the temperature at 100 kpc. We show the 
best-fitting model in Figure [S] and list the parameters in 
Table |21 The power-law exponent is a p = 0.00 ± 0.02, 
indicating an isothermal profile. This result is also 
consistent with that obtained for the Chandra data by 
iBuote fc Lewis! ()2004f) : the normalization of the temper- 
ature profile obtained from the XMM data is also within 
1.2cr of that obtained from Chandra. 

Unlike the low-quality Chandra temperature profile, 
the XMM temperature data are not formally consistent 
with the isothermal / power-law model. A better fit 
can be achieved by adding another degree of freedom 
to account for the curvature of the profile in the log-log 
plot. We experimented with several models and adopted 
a model consisting of two power-laws joined by an expo- 
nential cut-off term, 



Radius (kpc) 

Fig. 5. — Projected temperature profile. Open circles represent 
our XMM observation. The models are emission-weighted projec- 
tions averaged over each annulus (see JSJ. However, for display 
purposes we do not bin the models in the plot. The solid line is 
the emission weighted T(r)* model. The dashed line is the emis- 
sion weighted power-law. In the ratio plot the triangles represent 
the power-law and squares are the T(r)* model. 



residuals of the innermost and outermost data points are 
affected substantially. 

We adopt this profile as our default in our subsequent 
analysis. In § 03 we assess how different choices of tem- 
perature profiles affect our mass measurements. There 
we also assess the importance of the last data point for 
determining the best model. 

7. MASS ANALYSIS 

7.1. Gravitating mass 

When hydrostatic equilibrium is a suitable approxi- 
mation for the hot ICM the gravitating mass (M grav ) 
enclosed within a certain radius r can be inferred from 
the density and temperature of the emitting plasma 
ljMathewslll978t iFabricant et alJll980|) : 



M grav (<r) 



kr 



Gfirrip 



rT 



d\np d In T 
dh\r d In r 



(5) 



T(r)* = Tie 



-(&> 



T 2 (l 



(4) 



where ks is the Boltzmann's constant, G is the constant 
of gravitation, /i is the mean atomic weight of the gas 
(taken to be 0.62), m p is the atomic mass unit. We adopt 
the cusp (5 model for the gas density and the T{r)* model 
for the temperature for our fiducial analysis. Using these 
parameterizations we evaluate the gravitating mass at 



,3/2 



rf/ 2 )/2] 2 / 3 (see lLewis et aTl 



where Ti—± t 2 — ^.j(f-J and A;, r p , oil and 7 are re- 
spectively the normalization, scale radius and slope for 
the power-laws and the exponent of the power-law func- 
tions in the exponentials. This model provides a better 
fit (Figure [5] and Table |3J, though only the fractional 



weigh ted radii, r = [(', uui 

2003) within the annuli used for spectral analysis ($3} 
The resulting mass data points are displayed in Figure|SJ 
The error bars on the data points are the standard devia- 
tions derived from the gas density and temperature pro- 
files obtained from each of the 20 Monte Carlo error simu- 
lations (Q. Because the mass data points are correlated, 
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TABLE 3 

Mass profile fit parameters 



Profile 


X 2 C ± cr x 2 /dof 


r s (kpc) 


c 


a r„ ir (Mpc) 


M„„.(1O 14 M ) 


C200 


Af 2 oo(lO 14 M ) 


NFW 


12.7 ±7.3/5 


286 ± 27 


6.1 ±0.3 


1.74 ±0.05 


3.25 ±0.29 


4.6 ±0.3 


2.65 ±0.20 


N04 


2.4 ± 1.9/4 


206 ± 21 


7.0 ±0.3 


0.40 ±0.05 1.44 ±0.08 


1.84 ±0.33 


5.5 ±0.3 


1.70 ±0.24 


Power-law 


0.2 ±0.2/1 






1.84 ±0.18 


0.12 ±0.01 






Note. — 


For the N04 model r s corresponds to r_ 


2 . For the power-law model 


a corresponds to 


ol v and the 


M v i r value is the 



normalization of the model at 100 kpc. The Power-law model is fitted only to the inner three points. 
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Fig. 6. — Gravitating mass profile obtained using the cusp f3 
model for the gas density and the T(r)* model for the temperature. 
The gravitating mass is evaluated at radii corresponding to the 
annuli used for spectral analysis We display the best-fitting 

NFW (solid line), N04 (dashed line) and power-law (only for the 
three inner data points; dotted line) models. 



it is inappropriate to use the x 2 null-hypothesis proba- 
bility assuming uncorrelated errors to assess goodness-of- 
fit. Instead for each mass model we provide the standard 
deviation of \ 2 obtained from fitting a particular mass 
model (e.g., NFW) to the mass profile obtained for each 
of the 20 Monte Carlo simulations. This allows one to 
assess the relative goodness-of-fit between different mass 
models. 

First we examin ed whether the NFW profile 
ijNavarro et alJ 1199^ could provide a good description 
of the mass data. The NFW density profile is given by: 



p(r) = 



Pc(z)6 c 



(r/r a )(l + r/r s 



(6) 



where p c (z) is the critical density of the universe at red- 
shift z, S c is a characteristic dimensionless density defined 
as: 



1 pv 



3p c (z) ln(l + c)-c/(l + c) : 



(7) 



where r s is the scale radius (the point where the log- 
arithmic slope reaches a value of -2), p v ; lr is the mean 
density of a sphere enclosed in r„j r = cr s (where c 
is the concentration pa rameter) that is 104.7p c (z) (see 
iBrvan fe Normal! Tl 998). By integrating equation we 
obtain the expression of the total mass enclosed within 
a radius r: 



M NFW (< r) = A NFW 



ln(l + r/r s ) — 



r/r s 



1 + r/ r s 



(8) 



where Anfw — ^Pc{z)S c r^. A distinctive feature of the 
NFW profile is that the logarithmic slope of the density 
profile asymptotes to -1 at small radius corresponding to 
a logarithmic slope in the mass of 2. 

We plot the best-fitting NFW model in Fig. © (solid 
line) and list the derived parameters in Table |21 This fit 
is very good in the sense that the fractional residuals are 
< 5% for the central radial bin and for r > 100 kpc. The 
largest deviation is observed for r ~ 50 kpc where the 
fractional residual is ~ 13%. It is interesting that this 
radius corresponds to the largest center shift in the X-ray 
isophotes (©. The values of c and M vir obtained for the 
NFW fit are very consistent with those obtained from re- 
cent observ ations of relaxed clusters with Chandra and 
XMM (e.g.. iLewis et al.ll2003l IPointecouteau et alJl2005t 
Vikhl inTn et alJl200(tl . f These results are also consistent 
with those we obtained previously from the low-quality 
Chandra AO-3 data.) These c and M vir values agree 
extremely well with th e mean relatio n pred icted by the 
ACDM simulations of iBullock et alJ l|200lD . where we 
have extended their toy model up to virial masses ap- 
propriate for A2589. 

To obtain an estimate of the slope of the inner density 
profile without reference to the NFW model we fitted 
the three inner mass data points with a power-law (see 
Eq.|21 thin dotted line in the plot). The value we obtain is 
1.84±0.18 for the mass which corresponds to — 0.84±0.18 
for the density. This value is consistent with the inner 
logarithmic slope given by an NFW profile within the 
la error. The inner slope of -0.8 is also consistent with 
a direct solution of the radial Jeans equation for DM 
assuming an isotropic vel ocity dispersion tensor and a 
CDM phase-space density ( H ansen fc Stade]ll2005[) . 

Finally, we also examined the Sersic-like profile 
(jNavarro et alJ 120041 hereafter N04) that was recently 
suggested as a better parametrization for CDM halos, 
particularly within the very inner regions (less then few 
percent of the virial radiu s) where it is typically shal- 
lower than NFW (see also iGraham fc et alJ feOOS). The 
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N04 density profile is defined as follows, 
p(r) = p_ 2 exp 









(a, 


exp 






a 



(9) 



where x — r/r_2, and the quantities flagged with -2 refer 
to the point where the profile reaches a logarithmic slope 
of -2 (r_2 is the analog of r s in the NFW profile). The a 
exponent determines the bend of the profile about r_2- 
The enclosed mass obtained by integrating this expres- 
sion is: 

, , . f (2\(aV /a /3 2 c 
M NM {< r) = A NQi - exp - - 7 -, -x c 
a \a J \2 J \a a 



where A 



N04 



47T / o_ 2 r!_ 2 and 7(7/, A) = V V 



(10) 



the lower incomplete gamma function. 

We display the best-fitting N04 profile in Figure rep- 
resented by a dashed line. The fractional residuals are 
smaller than the NFW profile at all radii except the cen- 
tral data point where the residual is about 50%. The 
smaller \ 2 value obtained for the N04 arises primarily 
from the better fit of the data point r ~ 50 kpc. The 
inferred value of a = 0.40 ± 0.05 is quite large and in- 
comp atible with the mean value of O.f 72±0.032 for CDM 
halos l)Navarro et al.ll2~004|) . Consequently, the result we 
obtain for the concentration and virial mass using the 
N04 model must be interpreted with caution; i.e., the 
large a implies a density profile that is shallower in the 
center and steeper at large radii than CDM. 

Using our measurement of the gravitating mass ob- 
tained from the NFW model and the gas mass obtained 
by integrating the cusp (3 model we compute the gas frac- 
tion as a function of radius. For the central data point we 
obtain f gas = O.Of 4 ± 0.004 while it rises to 0.083 ± 0.002 
for our outermost mass data p oint. This value i s con- 
sist ent with values obtain ed bv [Sanders on et al.l (2003) 
and lVikhlinin et aTl l|2005fl for clusters of similar temper- 
ature. If we extrapolate our fits to the virial radius we 
obtain a total gas fraction f gas = 0.168 ± 0.009, consis- 
tent with the value of the baryon fraction Q,h/Vt m ~ 0.17 
measured by WMAP (|Spergel et alJl2003|) . 

7.2. The Ratio of Gravitating Mass to cD Stellar Light 

The radial profile of the ratio of gravitating mass to 
stellar light in the cluster provides a useful diagnostic 
for the radius where dark matter dominates the clus- 
ter mass budget. We shall only consider the optical 
light of the cD galaxy NGC 7647 since only sparse in- 
formation in the literature exists for a small number of 
other member galaxies. Fortunately, the U-band sur- 
face bri ghtness profile of the cD gala xy has been ana- 
lyzed by Malumuth & Kirshncr ( 1985) out to a radius of 
~ 130 kpc. They find that a King model is a good fit 
to the optical data with a core radius, r c = 0.98 h^g kpc, 
and a total U-band luminosity Ly — 2.3 x 10 11 Ly.o 
within the region analyzed. They find that a de Vau- 
couleurs model with r e — 32.6 kpc fits the optical 
data nearly as well. Therefore, we r epresent Lv(r) by 
a Hernquist model l)Hernauistl 1 19901) with scale radius 
a = r e /1.8153 = 18.0 kpc and total luminosity quoted 
above, since the projected Hernquist profile is a good 
approximation to a de Vaucouleurs profile. 

In Fig. [7| we plot M grav (< r)/Ly(< r): the ratio of 
gravitating mass, computed using the NFW model, to 



the integrated V-band luminosity of the cD, using the 
Hernquist model, as a function of radius. The M grav f Ly 
profile increases from 9.0 ± 2.3Mq/Lv,© for the inner 
data point to 93.5 ± 2.3M Q /Lv,o at the radius repre- 
senting the outer extent of the optical data. If we ex- 
trapolate the Hernquist model of the cD out to the ex- 
tent of the X-ray measurements, Mnr av /Ly increases to 
433±9M Q /L v ,o (gray area in Fig.0. This latter value 
is an upper limit to the cluster mass to light ratio be- 
cause we have not considered additional contributions 
from other member galaxies. 

There is no evidence that the M grav /Ly profile flat- 
tens at small radiu s such as we have observed for A2029 
llLewis et al.ll2003|) and for a sample of early-type galax- 
ies l|Humphrev et al.l 12006^1 . Since the flattening in the 
profiles in other systems occurs for r < 10 kpc, yet our 
central bin has outer radius r = 32 kpc, we do not ex- 
pect to detect this effect. High quality Chandra data 
would allow M glav /Ly to be probed on the small scales 
necessary for an inte resting comparison to other sy stems. 

We mention that iMalumuth fc Kirs hncr ( 1985) used 
their fit of the King model and a measurement of the 
central velocity dispersion to estimate M glav /Ly at the 
very center; i.e., within the King core radius (< 1 kpc). 
They obtain M grav / Ly rs 15 in solar units, converting 
their value to our redshift and cosmology. This value is a 
factor of 1.7 larger than we have measured at r ~ 20 kpc, 
although the discrepancy is not highly significant (2.6c). 
For a single burst stellar population with age ranging 
from 10-13 Gyr and metallicity ranging from 0.5-2 solar 
the U-band stellar mass-to-light ratio is only expected to 
take values from 5-10 in solar units, with values at the 
upper en d assuming a Sa lpeter IMF. (We have used the 
results of Maraston 199^ from updated model grids pro- 
vided by the author. ) The gravitating mass we measure 
for r ~ 20 kpc is similar to that we have measured at 
that radius in the cluster A2029 and a sample of ellip - 
tical galaxies l|Lewis eTaIll2003l: iHumphrev etai] l20061 . 
If the difference in masses estimated from the X-ray and 
optical methods is real, it likely implies the stellar mass- 
to-light ratio in the cD varies with radius. 

7.3. Dark matter 

We desire to extract the radial profile of the DM from 
the gravitating matter by removing the contributions of 
the known luminous mass components. Since the X-ray 
data provide a direct measurement of the hot gas den- 
sity, it is straightforward to compute the radial profile 
of M grav — M gas . Consequently, as a first step we fit- 
ted single-component models to this "DM" profile anal- 
ogously to our procedure in the previous section. 

The fitted parameters for the NFW, N04, and power- 
law models are listed in Table 0] Since the hot gas con- 
tributes < 10% to the total mass over the radial range 
fitted, it is not surprising that the parameters we infer 
are very consistent with those obtained when fitting the 
total gravitating matter (Table 0J. However, we empha- 
size that r v ; r and Af v i r (and hence c) obtained in this 
way refer only to the DM component. To obtain the 
true values that include the mass of the hot gas we add 
-Mgas to Mdm, compute a new r v ; r , and iterate. For 



http: / /www-astro. physics. ox. ac.uk/ ~maraston/ 
dia's_Stellar_Population_Models . html 
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Fig. 7. — Ratio of the total gravitating mass to the optical light 
of the cD galaxy. The data in the gray area are estimated from 
an extrapolation of the Hernquis t profile approximation of the De 
Vaucoulcurs profile obtained by Malumuth & Kirshncr ( 1985) to 
map the V-band light profile out to ~ 130 kpc. 



the NFW model this yields r vir = 1.77 ± 0.04 Mpc, 
Mvir = 3.46 ± 0.26 x 10 14 M Q , c = 6.8 ± 0.4, and a total 
gas fraction, f g = 0.16 ± 0.01. It should be remembered 
that these re-scaled values depend on the extrapolation 
of our model for p g out to r v ; r . 

Next we attempted to better isolate the DM profile by 
removing the contribution of the mass from the stars in 
the cD. To model the cD stars we employed the Hernquist 
parameterization of the y-band luminosity described in 
the previous section, which we now refer to as the H90 
model. We fitted two-component models, DM + stars, to 
the radial profile of M grav — Af gas , where the total mass 
in the stars is specified by the stellar mass-to-light ratio, 

M^/Ly. 

The results of the fits are listed in Table^J If M*/Ly is 
allowed to be a free parameter in the NFW+H90 fit, then 
an unphysically small value M^/Ly ~ 0.15 (solar units) 
is obtained. The NFW model itself is all that is required 
to describe the mass profile; adding in a separate stellar 
component only degrades the fit. 

We desired to explore the importance of adiabatic 
contraction of the DM halo arising from the early 
condensation of baryons in the central galaxy (e.g., 
iBlumenthal et al]ll986UGnedin et aUl2004j) . We refer to 
the modified NFW fit as "NFW*AC+H90" , indicating 
that adiabatic contraction only applies to the DM com- 
ponent 5 . Because the fitted M*/Ly is so small, the AC 
model does not have a substantial effect when M*/Ly is 
allowed to be a free parameter. Consequently, we also 

5 The adiabatic contraction code we used 
was made publicly available by Oleg G nedin at: 
http: / /www. astronomy.ohio-state.edu/ ~ognedin/contra/ 



tried fixing M*/Ly — 9.0 (solar units) corresponding to 
the gravitating mass-to- light ratio obtained for our inner 
data point, as well as being a reasonable estimate from 
stellar population synthesis models, as discussed in the 
previous section. With this M*/L v the NFW+H90 and 
NFW*AC+H90 fits are quite poor within r ~ 50 kpc, 
and the effect of adiabatic contraction on the derived pa- 
rameters for the DM is roughly equal to that arising from 
adding the stellar component; e.g., c changes from 6.3 
(NFW) to 5.2 (NFW+H90) and from 5.2 (NFW+H90) to 
4.3 (NFW*AC+H90). The adiabatic contraction of the 
dark matter produces a more centrally concentrated halo 
that deviates even more strongly with the data within 
r ~ 50 kpc (see Figure (SJ. 

Since the N04 model fitted to the total gravitating 
matter in fc l7.ll underestimates the mass at the center, 
it is expected that, unlike the NFW profile, it will al- 
low for a substantial contribution from a central stel- 
lar mass component. Indeed, the N04+H90 fit yields a 
stellar mass-to-light ratio, M*/Ly = 4.8 ± 1.9, which is 
consistent with expectations of single-burst stellar popu- 
lations synthesis models. Applying adiabatic contraction 
(i.e., N04*AC+H90) yields a smaller, though still reason- 
able, value, M*/Ly = 3.1 ± 0.7. However, the values of 
a ~ 0.5 obtained for these models, even larger than ob- 
tained for the fit to the gravitating matter, are even more 
inconsistent with t he values of ~ 0.20 required to repre- 
sent CDM halos ijNavarro et al J 1200^1 . Also, exactly as 
for the total gravitating matter, it is the data point near 
r ~ 50 kpc that is largely responsible for the better fit of 
the high-a N04 model with respect to NFW (see Figure 

8. ASSESSMENT OF SYSTEMATIC ERRORS IN THE MASS 
PROFILE DERIVATION 

In this section we provide estimates of the magni- 
tude of systematic errors on the parameters parame- 
ters we have deduced for the dark matter. We illus- 
trate the effects of systematic errors on two fiducial two- 
component models consisting of a stellar component and 
adiabatically compressed dark matter: NFW*AC+H90 
and N04*AC+H90. The results are summarized in Ta- 
bleEl The statistical error of the default model (A stat ) is 
also listed in the table. For each variant model involved 
in the systematic checks, we also obtain a corresponding 
statistical error (A sys ). If A sys >= A stat we also quote 
its magnitude with the associated best-fitting parame- 
ter shift in the table. For example, the concentration 
parameter for the NFW*AC+H90 model has a shift of 
0.4±0.4 in the A bkgmod column; i.e., A sys = A sta t in this 
case. We emphasize that none of the systematic checks 
that we explored can increase the inferred M^/Ly for 
the NFW*AC+H90 model to values 5-10 expected from 
the stellar population. 

In the following we describe the different checks per- 
formed: 

1. PSF: As stated in <2]we defined our annuli to have 
at least a width of 1'; in particular the central ra- 
dial bin is a circle of radius 30". With this choice 
the central circle encloses about 90% of the en- 
ergy of a point source. Furthermore, since the ra- 
dial temperature profile is observed to be nearly 
isothermal we expect that our measured "spectro- 
scopic" temperature should differ negligibly from 



DM radial profile of A2589 



11 



TABLE 4 
Model Fits to M„ r . 



- M a . 



Single-Component Models: Dark Matter Only 



Model 


M*/L v 


Xc ± CT X 2 /do} 


r-s(kpc) 


c 


a 


r„ ir (Mpc) 


M„ ir .(lO 14 M ) 


NFW 




12.4 ±7.1/5 


264 ± 26 


6.3 ±0.4 




1.65 ±0.05 


2.78 ±0.25 


N04 




2.4 ± 1.9/4 


193 ± 19 


7.1 ±0.3 


0.41 ±0.07 


1.37 ±0.08 


1.58 ±0.28 


Power-Law 




0.2 ±0.2/1 






1.82 ±0.18 




0.11 ±0.01 








Two-Component Models: Dark Matter 


+ cD Stars 






Model 


M*/L v 


Xc ± CT X 2 /dof 


r 3 (kpc) 


c 


Q 


r„ ir (Mpc) 


M mr (10 14 M Q ) 


NFW+H90 


(9) 


42.0 ± 18.0/5 


338 ± 21 


5.2 ±0.3 




1.74 ±0.05 


3.26 ±0.29 


NFW+H90 


0.15 ±0.02 


12.7 ±7.1/5 


263 ± 16 


6.2 ±0.3 




1.65 ±0.05 


2.78 ±0.25 


NFW*AC+H90 


(9) 


76.9 ±25.2/5 


428 ± 32 


4.3 ±0.3 




1.84 ±0.07 


3.84 ± 0.45 


NFW*AC±H90 


0.055 ± 0.001 


12.6 ±7.2/4 


264 ± 18 


6.2 ±0.4 




1.65 ±0.05 


2.78 ±0.25 


N04+H90 


(9) 


2.0 ±4.2/4 


190 ± 6 


6.6 ±0.1 


0.61 ±0.05 


1.25 ±0.03 


1.20 ±0.09 


N04+H90 


4.8 ± 1.9 


0.7 ±0.7/3 


190 ± 9 


6.8 ±0.2 


0.51 ±0.05 


1.30 ±0.04 


1.34 ±0.12 


N04*AC+H90 


(9) 


8.2 ± 7.1/4 


195 ± 6 


6.1 ±0.1 


0.81 ±0.08 


1.20 ±0.03 


1.04 ±0.07 


N04*AC+H90 


3.1 ±0.7 


0.8 ± 1.1/3 


191 ± 7 


6.8 ±0.1 


0.53 ±0.06 


1.29 ±0.04 


1.32 ±0.13 



Note. — For the power-law model M^ir is the value of the normalization at 100 kpc. Fixed parameters a re enclosed in 
parent hesis. "*AC" indicates the dark matter component is compressed adiabatically following the prescription of lGnedin et all 
120041) . Note that the virial radii and mass refer only to the dark matter component. 



Radius (arcmin) 



Radius (arcmin) 




50 100 
Radius (kpc) 



100 

Radius (kpc) 



500 



Fig. 8. — Models fitted to M gra v — M gas . (left panel) The solid line represents the pure NFW profile, the dashed line adds the 
stellar component (NFW±H90) with M*/Ly = 9.0 in solar units, and the dotted line applies adiabatic contraction to the NFW profile 
(NFW*AC±H90) with M*/L v = 9.0. (right panel) Same as left panel except the NFW profile is replaced with the N04 profile. 



the e mission-weighted tem perature within the bin 
fe.g.. Mazzott a et all2004l) . The consistency of our 
results with those obtained from the lower S/N, 
but higher resolution, Chandra data and 03, 
and with XMM and Chandra observations of other 
clusters (see below in i ilOjl , suggests that the larger 
XMM PSF has not biased our derived density and 
temperature parameters by amounts more than the 



estimated statistical errors. 

Background modeling {kbkgmod m Tabic OjJ: the 
major sources of uncertainty in the background 
modeling are the shape of the broken power-law 
that takes into account the instrumental back- 
ground at high energies and the level of the CXB 
component. We assess the influence of these contri- 
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TABLE 5 

Estimation of the systematic errors 



Parameter 


Bestfit 




^bkgmod 


NFW*AC±H90 

A-bndwidt A mod (T/p) 


^■exclpts (1 


6/2-7) 


Acen 


^■deproj 


c 


6.2 


±0.4 


±0.4 ±0.4 


-0.7 


±0.6 ±0.6/ - 1.0 


-0.6/ + 


0.6 




-2.4 




2.78 


±0.25 




+0.71 


-0.30/ ± 1.22 


±0.92/ - 


0.36 




+4.11 


M t /Ly 


0.055 


±0.001 


-0.054 


-0.051 


-0.044/ - 0.041 


-0.045/ - 


0.043 


-0.044 


-0.025 



N04*AC±H90 



Parameter 


Bestfit 




^bkgmod 


^■bndwidt A mod (T/p) 


A ea;c;p t s (l-6/2-7) 


Acen 


^■deproj 


c 


6.8 


±0.1 


±0.5 


±0.6/ -0.4 


±0.3/ + 0.5 


-0.3 


-1.0 


lyfyir 


1.32 


±0.13 


-0.21 


+0.16/ -0.10 


-0.50/ + 0.16 ±0.38 






a 


0.53 


±0.06 


±0.07 ±0.05 


-0.14/ + 0.34 


±0.26/ - 0.13 




+0.33 


M*/L v 


3.1 


±0.7 


-2.4 


-0.8 + 1.1 (-3.1, ±4.0)/ ±2.5 


+1.5 + 1.2/ - 3.1 


±2.1 


+3.3 + 2.2 



Note. — A B t a t is the ltr statistical error estimate. M v ir is expressed in units of 10 14 Mq. 



butions by a variation of ±20% in the slope of the 
high energy part of the broken power-law and in 
the normalization of the CXB component. We also 
set the parameters of the broken power-law to the 
values obtained fitting the spectrum of the out of 
field of view events. This last check almost always 
produced the major source of uncertainty, espe- 
cially because it caused the outermost data point to 
vary appreciably. However, this check does not pro- 
duce results that differ appreciably from the best 
fit value except for the mass to light ratio inferred 
for the N04*AC+H90 profile. 

3. We varied the band- width in which the spectrum in 
each annulus has been fitted (Abndwidt in TableEJ. 
We choose to restrict our analysis to the 0.5-5 keV 
and 1-5 keV bands. Only the NFW profile is af- 
fected by this check. 

4. We tried different models for the temperature and 
density profiles (that are shown separately in A mo d 
in Table |SJ) . The models we used are three for the 
temperature (the first three presented below) and 
one for the density. They are: 

(a) the power-law presented in Eq. [3J 

(b) a lognormal profile; 

(c) a profile that joins smoothly two power-laws: 
T(r) = (Pf 1 +pP 2 )- 1 / a , where Pi and P 2 are 
the power-laws (Eq. |3J) ; 

(d) a double (3 model with the /3's tied. 

It is worth noting that for both the DM parame- 
terizations the lognormal profile gives results very 
similar to the best fit. For the N04*AC+H90 
paramctrization the third profile does not allow for 
the presence of any central stellar component and 
the power-law temperature model predict a very 
high M/L. The different density parametrization 
has a noticeable impact on the NFW profile. 



5. We tried to exclude from the analysis either the 
first or the last annulus (A exc i pts in Table |SJ. We 
excluded them either only during the fitting of the 
dark matter profile or also in the temperature and 
density fitting. We show them separately in Table[5] 
as 1-6 (last point excluded) and 2-7 (first point ex- 
cluded). These points influence critically either the 
form of the dark matter profile (case 1-6) or the 
dominance of the central stellar component (case 
2-7). 

6. We repeated the analysis using annuli centered on 
the X-ray centroid of the cluster (A cent in Table|SJ). 
The parameters did not change much and their dif- 
ferences from the best fit values are smaller than 
the statistical errors. 

7. We analyzed the effect of the deprojection on the 
data {Adeproj in Tablc|SI) usin g the "onion-peeling " 
technique (e.g. iBuoteTfeoOOt iBuote et all l2003af) . 

The quality of the current data does not allow good 
constraints with this technique. 

9. NON-THERMAL PRESSURE SUPPORT 

Our analysis of the mass profile of A2589 inferred from 
the XMM data indicates that the NFW model provides 
a good description of the total gravitating mass profile. 
However, it is widely thought that the central dark mat- 
ter profile will be compressed adiabatically because the 
stellar baryons would have collapsed at ear ly times (e.g., 
IBlumenthal et alJll986llGnedin et alJl2004j) . As we have 
shown in H7.3l the adiabatic contraction model assuming 
a reasonable stellar mass-to-light ratio for the cD galaxy 
exceeds the mass data in the central r < 50 kpc. 

We investigate whether plausible additional pressure 
support from non-thermal processes could reconcile the 
XMM mass profile with the adiabatic contraction sce- 
nario. Since the most important non-thermal pressure 
should arise from turbulent motions and magnetic fields 
in the hot ICM, we shall focus our attention on them. 
We do not consider possible pressure support from cos- 
mic rays which, in principle, could be of the same magni- 
tude as that arising from magnetic fields. As pointed out 
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bv lLoeb fe Mad l)1994(l . the fact that S-Z measurements 
of clusters are not dominated by a synchrotron signal, 
the pressure from cosmic ray electrons should be much 
less than the thermal pressure in clusters. However, the 
pressure from cosmic ray pr otons may sti l l be i mportant. 

For simplicity we follow lLoeb fc Maol l) 19941) and as- 
sume, 

-^nonthcrm ^-^thcrm? (H) 

where a is a constant. This expression cannot be strictly 
valid since we expect little contribution from non-thermal 
pressure for r > 100 kpc. For our purposes we require 
only that da/dr can be neglected at the radii of inter- 
est. If we also assume that the equation of hydrostatic 
equilibrium remains a good approximation but with the 
thermal gas pressure, Pthorm, replaced by the total pres- 
sure, Pnontherm + -Pthernn then we have at any radius r, 



M(< r)NFW*AC+H90 ^ 

M (< r) NFW ' 



(12) 



where M(< t)nfw is the mass determined from fitting 
the NFW model to the total gravitating matter (Ta- 
ble Ej, and M{< r)NFW*AC+H90 is the mass obtained 
from the adiabatically contracted NFW model assuming 
M*/Ly — 9.0 in solar units for the cD galaxy (Table 0J|. 

For our inner mass data point (r = 20.2 kpc), where 
the discrepancy is largest (see Figure EJ, we obtain 
a = 0.93 using our best-fitting models. Note that the 
NFW*AC+H90 value is « 4a larger than the mass 
data point. For the second mass data point, r = 
45.3 kpc, we obtain a smaller value, a — 0.27, though 
the NFW*AC+H90 mass value remains ss 4cr larger than 
the mass data point. (Note that the NFW mass value is 
only 1. 5a larger than the second mass data point.) 

9.1. Turbulence 

If we assume that the only source of non-thermal pres- 
sure arises from random (isotropic) turbulent motions, 
-Pturb = 3Pc/( u turb)> tnen equations l(TT|) and ltT2*|l imply, 



Wturb 



'a 



k B T y 
3.26 keV/ 



1233 km s" 



(13) 



where the gas temperature corresponding to the inner 
radial bin (Tabled has been used. For our inner mass 
data point (r — 20.2 kpc) this gives, i>t U rb ~ l-3c s , where 
c s = 919 km/s is the adiabatic thermal sound speed of 
the gas. The generally undisturbed appearance of this 
cluster (SJ3J) argues strongly against a turbulent velocity 
this large. For the second mass data point, r = 45.3 kpc, 
we also have a very large velocity, Wturb ~ 0.7c s , where 
c s = 945 km/s corresponding to k B T = 3.45 keV. Cos- 
mological simulations generally find much smaller levels 
of turbulence in cluster cores, ftnrh = 0-1 — 0.3 c, s (e.g., 
iNagai et al]l2003t iFaltenbacher et al.ll2005t) In fac t, for 
lower mass clusters like A2589. lDolag et alJ ((2005) con- 
clude that only ps 5% of the total pressure can be ascribed 
to turbulent motions. 

A further check on the viability of turbulent pressure 
support is to verify whether, I, the eddy size correspond- 
ing to Kturb, is smaller than the length scale under con- 
sideration. In an approximate steady state the energy 
produced by turbulent motions must be dissipated on 
appropriate viscous length scales or else the gas will be 



rapidly heated. After setting the rate of turbulent energy 
generation equal to the rate of radiation energy loss, we 
solve for the eddy size, 



Pg 



k B T \ 15 

3.26kevJ V 2 - 6 x 10- 26 gcm- 3 
-l 

3.2 Mpc, (14) 



A(T, Z) 



1.7 x 10 -23 erg cm 3 s _1 

where p g is the gas density at r = 20.2 kpc and A(T, Z) is 
the plasma emissivity evaluated for the conditions of the 
inner annulus (Table ^ . The eddy size is vastly larger 
than the r ~ 20 kpc scale of the central region under 
consideration - and is even larger than the virial radius 
of ~ 1.7 Mpc. We conclude that pressure arising from 
random turbulent motions in the core cannot reconcile 
our measurement with the NFW*AC+H90 model. 

9.2. Magnetic Fields 

Now taking magnetic fields to supply all of the nec- 
essary non-thermal pressure (B 2 /8ir), we obtain from 
equations (|TD) and , 



B = ^fa 



3.26 keV 2.6 x 10" 26 gcm 



Pg 



0.5 



58 fiG, 



( 15 ) 

where, as above, we have taken T and p g appropriate for 
our inner mass data point, r = 20.2 kpc. For a = 0.93 
the required magnetic field for the inner mass data point 
is, B = 56/iG. Fields of similar magnitude have been 
suggested previously to explain X-ray mass discrepan- 
cies using ROSAT and Einstein data within the central 
< 1 kpc of the disturbed e lliptical galaxy NGC 4636 
l|Brighenti fc Mathewslll997fl . A field strength of 53 /xG 
was also suggested previously bv lLoeb fc Mad l|1994ft to 
explain the discrepancy between X-ray and lensing mass 
estimates in the core of the cluster A2218. 

But the typical magnetic fields in galaxy clusters range 
from 1 — 10/iG, much lower than requi red for interesting 
press ure support (e.g., for a review see lGovoni fc Ferettil 
2004). Larger fields have been measured in some clusters 
with strong radio sources. A particularly interesting case 
is A2029, which, though more massive, appears to be 
very relaxe d both within the core and on larger scales 
like A2589 l|Lewis et aUl2002fl . Unlike the radio-quiet 
A2589, the cluster A202 9 possesses a Wide- A ngle- Tail 
radio source from which lEilek fc Owenl l)2002|) estimate 
a field of 16/iG within 8 kpc of the center of A2029. This 
field does not provide interesting pressure support, and 
the strength likely d eclines rapidly with increasing radius 
l|Govoni et al]l2001ft . 

Although we recognize that there remain outstanding 
issues in the determinations of magnetic fields in clusters, 
current evidence clearly does not favor field strengths of 
« 56/iG in galaxy clusters. Therefore, we believe that 
magnetic pressure support is an unlikely explanation to 
reconcile the NFW*AC+H90 model with the XMM mass 
data we have presented for A2589. 

10. CONCLUSIONS 

Using a new XMM observation we have presented an 
analysis of the radial mass profile inferred from the prop- 
erties of the hot ICM of the radio-quiet galaxy clus- 
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ter A2589. We have confirmed the highly regular X- 
ray image morphology indicated by previous X-ray ob- 
servations possessing lower resolution an d/or lower S/N 
Ipavid et al.lll996HBuote fc Lewisll2004[) . The only no- 
table deviation of the X-ray image from concentric el- 
lipses is a si 16 kpc shift of the X-ray center of the 
R — 45 — 60 kpc annulus. The radial temperature pro- 
file is nearly isothermal unlike m ost X-ray regular clus- 
ters w hich display cool cores fe.g-. IDe Grandi fc Molendll 
2002). However, the metallicity does peak at a value near 
solar at the center and falls off with incr easing radius sim- 
ilar to that in cool-core clu sters fe.g.. iDe Grandi et al.1 
l200l iBohringer et all 120041. Overall, given the highly 
regular nature of the X-ray image and spectral proper- 
ties, and the fact it is bright, nearby, and possess no 
significant radio emission from the central galaxy, A2589 
is an especially good target for X-ray studies of its mass 
distribution. 

We find that the NFW profile fits the total gravi- 
tating matter (M grav ) very well over the region stud- 
ied (r — 0.015 — 0.25 r v i r , fractional residuals < 10%) 
with Cv ir = 6.1 ± 0.3 and M vir = 3.3 ± 0.3 x 1O 14 M 
(jvir = 1-74 ± 0.05 M pc) in excellent agre ement with 
the ACDM prediction (jBullock et al]l2001[h However, 
if we attempt to add a component to the mass profile 
representing the stellar mass of the cD galaxy with a 
reasonable stellar mass-to-light ratio, the fit is degraded 
substantially in the central ~ 50 kpc. Modifying the dark 
matter halo as the result of adiabatic contraction arising 
from the early condensation of stellar baryons in the cD 
(e.g. JBlumenthal et alll986HGnedin et all2004|) further 
degrades the fit. 

If instead we use the Sersic-like profile proposed by 
INavarro et alJ lj2004j) to represent CDM halos, then siz- 
able stellar mass-to-light ratios are implied that are rea- 
sonably consistent with predictions from single burst stel- 
lar population models. However, the inverse Sersic index, 
a ~ 0.5, obtained from the fits is a factor of ~ 3 higher 
than predicted; the dark matter profile in the core is in- 
ferred to be shallower than CDM. Hence, we are unable 
to obtain a good fit with a model consisting of CDM 
halo and a separate contribution from the stars in the 



cD galaxy with a reasonable stellar mass-to-light ratio. 

The good fit of the NFW profile to the total gravi- 
tating matter of regular galaxy clusters appears to be 
a common feature of X-ray studies. The bright, highly 
relaxed cluster A2029 follows the NFW profile all the 



way down to w 0.001r v 



d .1 12003ft . Dedicated 



studies of small samples of other mostly relaxed clus- 
ters with Chandra an d XMM l|Pointecouteau et alJl2005t 
Vikhl irdn et al]l2006fl also do not report significant devi- 
ations from the NFW profile arising from central stellar 
mass, except for some lower mass, group-sca le objects 
(|Vikhhnin et alJl2006t iGastaldello et al]l2006ft . 

Since it seems unlikely that deviations from hydrostatic 
equilibrium in the hot ICM have conspired in the same 
way in all of these (regular) clusters to produce a mass 
profile consistent with NFW in the center, we conclude 
that the adiabatic contraction scenario does not appear 
to describe the formation of X-ray clusters. In partic- 
ular, for A2589 we have estimated the amount of non- 
thermal pressure support in the hot ICM from random 
turbulent motions and magnetic fields and conclude that 
neither source, especially turbulence, can reconcile the 
XMM mass profile with the adiabatic contraction sce- 
nario assuming a reasonable stellar mass-to-light ratio 
in the cD. We suggest that X-ray observations of A2589 
and other relaxed clusters favor the scenario where pro- 
cesses during halo formation, such as the heating of the 
dark matter by dynamical friction with member galaxies, 
counteracts adiabatic compression and leads to a total 
gravitating mass profile consist ent with the pure NFW 
profil e fe.g.. lEl-Zant et afll2004l see also lLoeb fc Peebles! 
2003). 
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